Lifetime of dynamic heterogeneities in a binary Lennard-Jones mixture 
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A four-time correlation function was calculated using a computer simulation of a binary Lennard- 
Jones mixture. The information content of the four-time correlation function is similar to that of 
four-time correlation functions measured in NMR experiments. The correlation function selects a 
sub-ensemble and analyzes its dynamics after some waiting time. The lifetime of the sub-ensemble 
selected by the four-time correlation function is calculated, and compared to the lifetimes of slow 
sub-ensembles selected using two different definitions of mobility, and to the a relaxation time. 

PACS numbers: 61.43.Fs,64.70.Pf 



The origin of the non-exponential relaxation found in 
supercooled liquids has been studied extensively in the 
last ten years. Two possibilities exist P, Q ■ Either all 
the particles undergo non-exponential relaxation (homo- 
geneous scenario) , or the relaxation of each particle is ex- 
ponential and there is a large variation in the relaxation 
time of the particles (heterogeneous scenario). There 
have been many simulations 

iliilliillflia and 

experiments ^2 Q> HM Il6l| which imply heteroge- 
neous relaxation. The heterogeneous relaxation scenario 
suggests that the particles in a supercooled liquid can be 
categorized by their relaxation time. The particles with 
the shortest relaxation times are referred to as "fast" par- 
ticles, and the particles with the longest relaxation times 
are "slow" particles. One important question is the life- 
time of the dynamic heterogeneities, i.e. how long does 
a fast particle remain fast and a slow particle remain 
slow? The first part of this question was considered in 
one of the early simulational investigations of dynamics 
heterogeneities : the lifetime of fast particles has been 
found to be much shorter than the a relaxation time. It 
should be noted that experiments are usually sensitive to 
slow particles and thus simulational investigation of the 
slow particles lifetime is also important; however, to the 
best of our knowledge, lifetime of slow particles has been 
studied only in two dimensions where it has been found 
to be comparable to the a relaxation time ^3 • Here we 
study the lifetime of slow particles using an approach in- 
spired by one of the experimental protocols. Our study 
is complementary to recent investigations of the spatial 
correlations of the slow particles [lO, ■ 

The lifetime of dynamic heterogeneities has been mea- 
sured in a reduced four-dimensional nuclear magnetic 
resonance (NMR) experiment by monitoring parts of a 
four-time correlation function. The general idea of the 
experiment has been lucidly explained by Heuer |l8j |: 
one can define a filtering function /(ti,t2) such that 
(/(iij ^2)) selects particles which are slow over a time in- 
terval Ati2 = t2 — ti. Thus, {f{ti,t2)f{t3,t4)) selects 
particles which are slow over time intervals Atu and 
Ai34 = ti — t^. The two time intervals are separated 
by a waiting time = — t2- For small t^, the re- 
laxation of the slow sub-ensemble remains slow, but for 
large enough t^, the relaxation of the slow sub-ensemble 



is the same as the relaxation of the full ensemble. The 
lifetime of the slow ensemble is related to the minimum 
tu, such that the average relaxation time of the slow sub- 
ensemble returns to the average relaxation time of the full 
ensemble. Bohmer et al. j3| used this idea to investigate 
Ortho-Terphenyl (OTP) at 10 K above Tg ^ 243 K. Us- 
ing a pulse sequence they selected a set of particles which 
did not rotate appreciably over a time interval Aii2, i.e. 
a slow sub-ensemble. The particles were then allowed 
to evolve during a time interval t^. Finally they mea- 
sured what fraction of the slow sub-ensemble were still 
slow over a time interval At^^. The characteristic time 
for the slow sub-ensemble to remain slow was found to be 
comparable to the average relaxation time of the full en- 
semble. This is in a stark contrast with results obtained 
for OTP by Ediger's group[ll[l|: at Tg + 4 K the life- 
time of the dynamic heterogeneities was found to be 6 
times longer than the a relaxation time and at Tg -I- 1 K 
it was 100 times longer! Ediger's findings could, however, 
be compatible with the NMR result if strong temperature 
dependence of the lifetime sets in close to Tg . 

The procedure used in this work to measure the life- 
time of dynamic heterogeneities is closely related to the 
NMR approach described above. We use a four-time 
correlation function to select a slow sub-ensemble, and 
monitor the relaxation and the lifetime of the slow sub- 
ensemble. The four-time correlation function selects a 
sub-ensemble without any explicit definition of mobility, 
thus it is not clear which particles are contributing to 
the four-time correlation function. To identify these par- 
ticles we use different definitions of mobility to select sub- 
ensembles whose relaxation is similar to the sub-ensemble 
selected by the four-time correlation function. Finally, we 
measure the lifetime of these slow sub-ensembles. 

To investigate the lifetime of dynamic heterogeneities 
we use the trajectories generated by an extensive Brow- 
nian Dynamics simulation study of a 80:20 mixture of 
a binary Lennard-Jones fluid [l3- Briefly, the potential 
is given by Vap = de^^ [(ca/j/r)^^ - (aap/r)^], where 
a,/3 e {A,B}, and eAA = 1-0, e^s = 1-5, esB = 0.5, 
aAA = 1-0, (TAB = 0.8, and aBs — 0.88. A total of 
N — Na + Nb — 1000 particles were simulated with a 
fixed cubic box length of 'd.AaAA- AH the results are pre- 
sented in reduced units where uaa and eaa are the units 



2 



of length and energy, respectively. The system was sim- 
ulated at temperatures T = 0.44, 0.45, 0.47, 0.5, 0.55, 
0.6, 0.8 and 1.0. A long equilibration run, and two to 
eight production runs were performed at each tempera- 
ture. The equilibration run was at least as long as the 
production runs. The presented results are the average 
of the production runs. The characteristics of this glass- 
forming liquid has been extensively studied TO^ 20, . 
The details and the results of the Brownian dynamics 
simulation are given elsewhere |2^ . In particular, we 
found that a relaxation times. Fig. |21 follow a power- 
law temperature dependence in the temperature range 
0.47 < T < 0.8 and deviate from this power-law depen- 
dence for T < 0.47. This is similar to earlier findings 
using Newtonian 1^, '2T| and stochastic dynamics [2Cj. 

To examine lifetime of dynamic heterogeneities we fol- 
low the procedure discussed above: we use a filtering 
function f{ti,t2) = e*i-('"j(*2)-'"^(*i», where rj{t) is the 
position of particle j at time t. Thus (/(ii, ^2)) is the in- 
coherent intermediate scattering function Fs{q;t2 — ti). 
For all the calculations, q is set to a value around the 
first peak in the AA {q = 7.25) or BB [q = 5.75) par- 
tial structure factor for Jv[^ and respectively. The 
four-time correlation function is defined as follows: 

{f{tiM)f{hM)) 

J_ ^ giq-(r3(t2)-r3(ti))giq-(rj(t4)-rj(t3)) \ 




(1) 



p'q-(rj(t2)-rj(ii)) 



where a G {A, _B}. The normalization of the correlation 
fimction is such that if = t^, then A^" = 1.0. For small 
tw = ts — t2, the relaxation of the slow sub-ensemble 
remains slow, but for large enough t^^ the relaxation of 
the slow sub-ensemble is the same as the relaxation of 
the full ensemble. 

We fix the first time interval, Aii2 = ^2 — ^i, to be 
equal to 3tq where is the a relaxation time (tq is 
defined by the usual relation Fs{q,Ta) = e~^). This is 
comparable to the longest time intervals Ati2 used to 
select a slow sub-ensemble in the NMR experiment of 
Bohmer et al. Note that the time Ati2 = is well 
past the plateau region of the mean squared displace- 
ment, and is longer than what has been used in previ- 
ous simulational investigations which examined dynamic 
heterogeneities |^^|^. The second time interval, the 
waiting time tw = ^3 — ^2, is varied. Finally, for a given 

M°'{q,tw,t) EE 7W"(q,0,3T„,3T„ -f i^,i-t-3T„ +t,„) 

is calculated as a function of time t {i.e. as a function 
of the last time interval, At34 = ~ t^). AI^{q,tw,t) 
is shown in Fig. ^ for several waiting times. Notice that 
if tw = 0, then M''{q,tw,t) = (g,3Ta +t)/i^f (g,3T„). 
Also, M°'{q,tw,t) converges to F"{q,t) as the waiting 
time increases. The lifetime of the sub-ensemble mea- 
sures how long it takes for this convergence to occur. 
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FIG. 1: F^{q,t) (solid line) and M^{q,t^,t) (dashed lines) 
for = 0, 5, 50, 250, 500, and 1000 at T = 0.45 listed in order 
from the longest relaxation time to the shortest relaxation 
time. (Insert) H^{q,t^,t)/H:^^^{q,0) for = 0, 5, 50, 250, 
500, and 1000 at T = 0.45. 



We define the lifetime of dynamic heterogeneities 
as the waiting time for which the difference between 
M"{q,tw,t) and F^{q,t) is equal to of its value at 
short times. The exact procedure is as follows: As shown 
in Fig. n for t^ > there is an initial decay of M" to a 
plateau region, then M" decays to zero after the plateau. 
This is in contrast to the t^ = case where there is 
no initial decay to a plateau. Since we are interested 
in the relaxation after the plateau, M"(5, 0,i) is mul- 
tiplied by a temperature dependent factor C'(T) so that 
C{T)M°-{q, 0, ic) = ^r(g, tc) where tc is at the beginning 
of the plateau region of F". The choice of tc affects the 
results slightly, with a larger tc leading to a somewhat 
longer lifetime. However, the choice of tc does not affect 
any of the conclusions of this work. We calculate 



H''{q,tw,t) 



C{T)M^{q,tw,t)-F^{q,t) tw^O 
M-^{q,tw,t)-F^{q,t) t„ > 



(2) 

and determine the lifetime as the waiting time when the 
peak value of 77" (g, t^, , t) , is a factor of e smaller than its 



til 



value, i.e. 



H^ax{q,T-f,)/H^^^{q,Q) = e 1, where 



^maxil^^w) is the maximum value of H°'{q,tw,t). 

Shown in Fig. |21 is the temperature dependence of the 
lifetime of the slow sub-ensemble selected by A^" and 
for comparison the a relaxation time. Notice that the 
lifetime is not longer than the a relaxation time. The life- 
time increases faster with decreasing temperature than 
the a relaxation time except at the lowest temperatures 
studied where it has the same temperature dependence 
as the a relaxation time. 

An advantage of a computer simulation is that the tra- 
jectories of individual particles can be followed through- 
out the simulation. This allows us to try to identify a slow 
sub-ensemble which is a major contribution to the four- 
time correlation function, i.e. the sub-ensemble selected 
by f{ti,t2). To this end we have defined the mobility 
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FIG. 2: The characteristic lifetime found using the four-time 
correlation function (a), by using Ui (Eq. to define the 
mobility (0), and by using Si (Eq. (0J) to define the mobility 
(o), compared to the Q-relaxation time (dashed line). 



(Ti(At) of a particle i over a time interval At as: 

a.{M)^\v.,{t)-v.,{h)\\ (3) 

where the bar denotes an average over time t £ (ti, ii + 
At) . A particle is defined as slow over a time interval 
At if ai (At) is less than a cutoff value r1^^ . These are the 
particles which stay closest to their position at ti during 
the whole time interval At. 

To make a connection with the four-time correlation 
function study we fix At — 3tq,. Next, the incoherent in- 
termediate scattering function, Ff^^^^q^t), is calculated 
for the slow particles after a waiting time t^ has elapsed. 
^siow shown in Fig.|31for different values of r^„t, and 
is compared to and M"^. Note that F^^^ and 
are calculated for the same waiting time tw = 0.2. For 
a large cutoff r^„( , the sub-ensemble behaves like the full 
ensemble. For smaller values of r^„j, the average relax- 
ation time of the slow particles is longer than the average 
relaxation time of the full ensemble. For a small enough 
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FIG. 3: F^{q,t) (dashed line), F^^^{q,tn,,t) (dotted lines) 
for rL* = 0.05, 0.03, 0.025, 0.02, 0.015, 0.014, 0.013 listed 
from left to right, and M^{q,t^,t) (solid line) for — 0.2. 



cutoff, F^i^^{q,t) « M°'{q,t^,t). The size of the cutoff 
needed to achieve this equality depends on tw and the 
time interval used to identify the slow particles. For the 
temperature shown in Fig. ^ F^^^^{q,t) « M^(g, 0.2,t) 
for r^„( = 0.015. This cutoff corresponds to the 0.075% 
slowest particles. As increases, the value of r^^^ re- 
sulting in F"ig^ « M" also increases. For the higher 
temperatures, it was not possible to find a value of r^^^ 
so that F°[^^ « Af" for short waiting times. 

The characteristic lifetime of the slow particles To- can 
be calculated using the algorithm described above (note 
that now we do not need the correction factor C{T)). The 
temperature dependence of the characteristic lifetime of 
the slow sub-ensemble is shown in Fig.|21 The cutoff was 
chosen so that on average the 10% slowest particles were 
used in the calculation. The choice of the cutoff has little 
effect on the lifetime, as long as a sub-ensemble with a 
relaxation time longer than the average relaxation time 
of the full ensemble is identified. The lifetime calculated 
by identifying the slow particles is always equal to the a 
relaxation time to within the uncertainty of the data. 

Refs. 0)13 used the following measure of the mobility 



5,(At) = |r,(t2)-r,(ti)| 



(4) 



where At = t2 — ti. We defined a slow sub-ensemble as 
the 10% with the smallest Si{3Ta), and calculated F"i^^ 
for this sub-ensemble. Again, the average relaxation time 
of the sub-ensemble was longer than the average relax- 
ation time of the full ensemble. The lifetime of the sub- 
ensemble defined using the second definition of the mo- 
bility, Tg, is equal to the a relaxation time to within the 
uncertainty of the data except for the A particles at the 
highest temperatures examined in this work (see Fig.[2J). 

To try to understand why both definitions give similar 
results, it is illustrative to examine the relaxation of dif- 
ferent subsets of particles chosen by the two definitions 
of mobility. Let S be the set of particles selected using ai 
as the definition of mobility, and V be the set of particles 
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FIG. 4: Fsiow for different different sub sets of particles, SPiT) 
(dashed dotted), S — V (dotted), T> — S (dashed), compared 
to (solid). See text for definition of sets S and T). 



selected using 5i as the definition of mobility. Figure 01 
compares for 5 n P, 5 - P, and 2? - 5 to for 

T = 0.55. The relaxation of the particles which are in set 
S but not T), or are in set T) but not <S, is similar to the 
relaxation of the full ensemble, but the particles which 
are in both sets have a longer relaxation time. Thus, the 
two definitions of mobility give similar results since they 
both are able to select the particles whose average relax- 
ation time is longer than the average relaxation time of 
the full ensemble. 

In conclusion, we used a four-time correlation function 
to select a slow sub-ensemble and analyze the dynamics 
of the slow sub-ensemble. The lifetime of the slow sub- 



ensemble selected by the four-time correlation function is 
not longer than the a relaxation time. On approaching 
the lifetime increases faster with decreasing temperature 
than the a relaxation time . Closer to Tc (beginning 
approximately at the temperature at which deviations 
from mode-coupling-like power laws appear) the lifetime 
follows the temperature dependence of the a relaxation 
time. We also identified two other slow sub-ensembles 
whose average relaxation time is longer than the average 
relaxation time of the full ensemble using two different 
definitions of mobility. The essential sub-ensemble, the 
sub-ensemble chosen such that F^^^^ ss M", consists of 
the particles which stay closest to their position at ti 
over the time interval At = t2 — ti, and are still close 
to their position at ti. This suggests that the slow sub- 
ensemble are the particles which are confined to their 
cage over the time interval At. The lifetime of the slow 
sub-ensemble depended on the definition of mobility. If 
ai was used to define mobility, the lifetime was equal to 
the a relaxation time at all temperatures. If 6i was used 
as the definition of mobility, the lifetime was equal to 
the a relaxation time except for the A particles at the 
highest temperature studied, in which case the lifetime 
was less than the a relaxation time. 

Our findings qualitatively agree with NMR results of 
Bohmer et al. [TJ. Note, however, that there is a signifi- 
cant difference in the temperature of simulations and ex- 
periments: the simulations have been performed slightly 
above Tc whereas the experiments were done well below 
Tc- Thus, direct comparison of the two sets of results 
is impossible. The same comment applies, however, to 
almost all simulational studies of glassy dynamics. 

This work was supported by NSF Grant CHE 0111152. 



[1] M. Ediger, Annu. Rev. Phys. Chem. 51, 99 (2000). 
[2] R. Richert, J. Phys. Cond. Matt. 14, R703 (2002). 
[3] W. Kob et ai, Phys. Rev. Lett. 79, 2827 (1997). 
[4] J. Qian and A. Heuer, Eur. Phys. J. B 18, 501 (2000). 
[5] J. Qian, R. Hentschke, and A. Heuer, J. Chem. Phys. 

110, 4514 (1999). 
[6] A. Heuer and K. Okun, J. Chem. Phys. 106, 6176 (1997). 
[7] C. Donati et al, Phys. Rev. Lett. 80, 2338 (1998). 
[8] B. Doliwa and A. Heuer, Phys. Rev. Lett. 80, 4915 

(1998) . 

[9] P.H. Poole, C. Donati, and S.C. Glotzer, Physica A 261, 
51 (1998); C. Donati et al, Phys. Rev. E 60, 3107 (1999). 
[10] N. Lacevic and S.C. Glotzer, J. Phys. Cond. Matt. 15, 
S2437 (2003). 

[11] N. Lacevic et al, J. Chem. Phys. 119, 7372 (2003). 
[12] R. Bohmer et al, Europhys. Lett. 36, 55 (1996); R. 

Bohmer et al, J. Chem. Phys. 108, 890 (1998). 
[13] M.T. Cicerone and M.D. Ediger, J. Chem. Phys. 103, 

5684 (1995). 

[14] C. Wang and M. D. Ediger, J. Phys. Chem. B 103, 4177 

(1999) . 

[15] F. R. Blackburn et al, J. Non-Cryst. Solids 172-174, 
256 (1994). 



[16] K. Schmidt-Rohr and H.W. Spiess, Phys. Rev. Lett. 66, 
3020 (1991). 

[17] D.N. Perera, P. HarroweU, J. Chem. Phys. Ill, 5441 

(1999); B. Doliwa, A. Heuer, J. Non-Cryst. Solids 307- 

310, 32 (2002) 
[18] A. Heuer, Phys. Rev. E 56, 730 (1997). 
[19] W. Kob and H.C. Andersen, Phys. Rev. E 51, 4626 

(1995); Phys. Rev. E 52, 4134 (1995). 
[20] T. Gleim, W. Kob, and K. Binder, Phys. Rev. Lett. 81, 

4404 (1998). 

[21] W. Kob, in Slow Relaxations and Nonequilibrium Dy- 
namics in Condensed Matter; J.-L. Barrat et al., eds. 
(EDP Sciences-Springer V erlag, 2003). 

[22] G. Szamel and E. Flenner, [cond-mat /0407537| (accepted 
for publication in Europhys. Lett.). 

[23] A similar definition was used below Tc by K. VoUmayr- 
Lee, W. Kob, K. Binder, and A. Zippelius, J. Chem. 
Phys. 116, 5158 (2002). 

[24] The a relaxation time and the lifetime of the slow sub- 
ensemble were both defined by the condition that an ap- 
propriately defined function decayed to an arbitrary value 
of e-^ If, instead of e'^ , 0.25, 0.35 or 0.5 were used, the 
same trends were observed. 



